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Abstract — A semi-analytical method evaluates the error prob- 
ability of DPSK signals with intrachannel four-wave-mixing 
(IFWM) in a highly dispersive fiber link with strong pulse 
overlap. Depending on initial pulse width, the mean nonlinear 
phase shift of the system can be from 1 to 2 rad for signal-to-noise 
ratio (SNR) penalty less than 1 dB. An approximated empirical 
formula, valid for penalty less than 2 dB, uses the variance of 
the differential phase of the ghost pulses to estimate the penalty. 
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I. Introduction 

DIFFERENTIAL phase-shift keying (DPSK) signal has 
been studied widely recently for long-haul lightwave 
transmission systems [l]-[4]. In additional to the 3-dB receiver 
sensitivity improvement to on-off keying (OOK), high-speed 
DPSK signal also has larger tolerance to fiber nonlinearities 
than OOK signal [5]. Most DPSK experiments use return-to- 
zero (RZ) short pulse and launch a constant-intensity pulse 
train with phase modulated to each RZ pulse. 

For 40-Gb/s signal in dispersive fiber, each RZ pulse broad- 
ens very fast by chromatic dispersion and overlaps with each 
other The pulse-to-pulse interaction gives intrachannel cross- 
phase modulation (IXPM) and four-wave-mixing (IFWM) 
[6], [7]. As a constant pulse train, IXPM induces identical 
phase modulation and timing jitter to all pulses and does not 
affect DPSK signal. However, IFWM adds ghost pulses to 
each DPSK RZ pulse [8]-[ll]. DPSK signal also has higher 
tolerance to IFWM than OOK signal [11]. 

When the IFWM induced ghost pulses are evaluated numer- 
ically, the error probability of DPSK signal can be calculated 
semi-analytically. This letter studies the statistical properties of 
IFWM in more detail using a method similar to [9], [11]. Both 
the error probability and signal-to-noise ratio (SNR) penalty 
are calculated. 

II. Statistics of IFWM 

If the signal launched to the fiber link is Gaussian pulse train 
with initial 1/e-pulse width of Tq or full-wide-half-maximum 
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(FWHM) pulse with of IMTq, for a bit-interval of T, the fcth 
pulse is Uk = ^fccxp [-(t- fcr)^/2To], where Ak = ±Aq 
is phase modulated by either or tt. From [8], [9], [11], the 
peak amplitude of the ghost pulses induced from IFWM, from 
the m, n, and {m + n)\h pulses to the pulse at t = is 
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where m ^ and n 7^ 0, 7 is the nonlinear coefficient and a 
is the attenuation coefficient of the fiber, /32 is the coefficient 
of group velocity dispersion, and L is the fiber length per span. 
Here, we exclude both IXPM with either m = or n = and 
self-phase modulation with m = n = 0. 

Figure [fla) shows the distribution of the normalized com- 
plex electric field of Auq/Aq with the unit of radian. Fig- 
ure [flb) shows the distribution of the peak phase shift of 
^{Auq}/Aq versus ^{Aui}/Ai between two consecutive 
time intervals, where Aug and Aui are the peak-amplitude 
of ghost pulses and 3{} is the imaginary part of a com- 
plex number The ghost pulses of Auq and Aui include all 
contributions of —8 < m,n,n2 + n < 8 for a 16-bit DPSK 
signal with about 64,000 combinations. Figure^a) is for ghost 
pulse at the center bit and Fig. [Qb) is for the center two 
bits. If the pulse amplitude of |j4o| is significantly larger than 
the IFWM ghost pulses, 3{Amo}Mo and 3{Aui}/Ai give 
approximately the phase shift [11]. 

Figures [Q are obtained for an A^-span fiber link with L = 
100 km of fiber per span with a normalized launched power of 
unity mean nonlinear phase shift of ($nl) — N-fLcsPo = 1 
rad, where Pq = y^jylo pTg/T is the launched power, and 
LcS ~ l/a is the effective nonlinear length per span. N 
identical fiber spans are repeated one after another with 100% 
dispersion compensation at the end of each fiber span. For 
arbitrary fiber link configuration, N instead of one integration 
of Q are required. IFWM ghost pulses add coherently span 
after span as the worst case. The fiber link has an attenuation 
coefficient of a = 0.2 dB/km. With bit interval of T = 25 
ps, DPSK signal has a data rate of 40 Gb/s. The initial pulse 
width is To = 5 ps, for a duty cycle of about 1/3. The fiber 
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Fig. 1. (a) The distribution of tiie complex electric field of Auo/Aq. (b) 
The distribution of 3{A«o}/^0 versus S{Ami}/Ai for two consecutive 
bit interval. 



dispersion is P2 = —22 ps^/km, corresponding to D = 17 
ps/km/nm at the wavelength of 1.55 /xm for standard single- 
mode fiber 

The distribution of Fig. ^a) is very irregular and has 
significant discrepancy with Gaussian distribution. Similar to 
Wei and Liu [11, Fig. 2], Figure [flb) is only symmetrical with 
respect to x ~ y. The difference of Fig.^b) with [11] is for 
a lossy instead of lossless fiber. With span by span dispersion 
compensation and for lossy fiber, 50% precompensation of 
dispersion increases instead of reduces IFWM. The phase of 
Q{Auq}/Aq is correlated with ^{Aui}/Ai with a correlation 
coefficient of about 0.58. For N repeated identical fiber spans. 
Figures ^ are valid for single- and multi-span systems with 
('I'nl) = 1 rad. Note that both Auq/Aq and Aui/Ai are zero 
mean. Not shown in Figs. [2 the real parts of di{Auo}/Ao and 
di{Aui}/Ai have a correlation coefficient of about -0.54. 

III. Error Probability for DPSK Signals 

The error probability of DPSK signals with IFWM is 
difficult to find analytically. From Figs. [2 the distribution of 
the IFWM induced ghost pulses is not Gaussian distributed. 
With the distribution of Fig. [2 the error probability of DPSK 
signal with IFWM can be calculated semi-analytically. 

Assumed for simplicity that the transmitted phases at t = 
and i = T are identical and, without loss of generality, the 
transmitted signals are Es{t) = Es{t — T) = Aq > 0. 
With optical amplifier noise of n{t), ignored the constant 
factor of interferometer loss and photodiode responsivity, the 
photocurrent is [12] 

i{t) = \2Ao + Aui + Auq + n{t) + n{t - T)\'^ 

-|Ami- Auo+n(t)-n(t-T)p. (2) 

A decision error occurs if i{t) < 0. 

Given Auq and Ami, the two terms in (|2j are independent 
of each other and have a noncentral chi-square distribution 
[13, pp. 41-44]. Each term of (|2} has the same noise variance 
of 4ct^ where i?{|n(i)p} — 2ct^ with a^ as the noise variance 
per dimension. The noncentralities of the two terms of (|5} are 
\2Af) + Aui + AuqP and |Ami — AwqP, respectively. From 
[13, App. B] [14], the probability of i{t) < is equal to 
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Fig. 2. The error probability of DPSK signals with IFWM ghost pulses as 
a function of SNR ps . 



where Q(-, •) is the Marcum Q function and 
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where ps = Aq/2<t'^^ is the SNR without taking into account 
the ghost pulses. Evaluated simi-analytically using the distri- 
bution of Figs. [2 the error probabihty is equal to 



Pe= E {pe{Auo, Aui)} , 



(6) 



where E{-} denotes expectation. 

When the sequence of Ak is changed to {—!)'' Ak with all 
odd positions changing sign, from Q, Auq remains the same 
but Aui changes sign. As Auo/Ao and Aui/Ai remain the 
same for Fig.fHb), the error probability for the case with Ao = 
—Ai is the same as that of (|6}. 

Figure |2] shows the error probability as a function of SNR 
Ps for DPSK signal with IFWM induced ghost pulses. The 



error probability without IFWM of 



[13, Sec. 5.2.8] is 



also shown for comparison. The semi-analytical formula of 
(|6j with (|3j is used to calculate the error probability based on 
IFWM ghost pulse distribution of Figs. ^ 

Figure|3]shows the SNR penalty forpe = 10^^ as a function 
of the mean nonlinear phase shift of (^nl)- In additional to 
the SNR penalty corresponding to Fig.|2j Figure |3l also shows 
the penalty when the initial pulse width is Tq = 7.53 ps for a 
duty cycle of 1/2. For 1-dB SNR penalty, the mean nonlinear 
phase shift must be less than 1.25 and 1.80 rad for initial pulse 
width of To = 5 and 7.53 ps, respectively. The SNR penalty 
is smaller for small initial pulse width of To = 5 ps. 

If the IFWM induced ghost pulses are assumed as Gaussian 
distributed, the noise increases to n{t) + Auo at t = 0. 
However, the SNR of \Ao\'^ /E{\n{t) + AuoP} cannot be 
used directly to find the error probability due the correlation 
between the ghost pulses of Auo and Aui . The phase variance 
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Fig. 3. The SNR penalty versus mean nonlinear phase shift of {'I'nl)- The 
curves are from empirical formula. 



of Aq + n{t) is equal to about l/2ps [13, Sec. 5.2.7]. The 
variance of the differential phase is 1/ps and should compare 
with 
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due to IFWM ghost pulses. The SNR penalty can be empir- 

where 20 is the 



20^2,) 



ically estimated as — 10 • logj^Q (l 

SNR for a DPSK error probability of 10^^. This empirical 
formula finds the SNR penalty based on the variance of the 
differential phase. Figure |3] also shows the SNR penalty from 
the empirical approximation. For SNR penalty less than 2 dB, 
the empirical approximation underestimates the SNR penalty 
by less than 0.25 dB. For SNR penalty larger than 2 dB, the 
approximation overestimates the SNR penalty. 

The empirical approximation still requires extensive nu- 
merical calculations to find the variance of a^g from many 
combinations of bit sequence. The semi-analytical method 
needs one further step to find the error probability of (|3j for 
each term of Amq and Ami, and then averaging of (|6jl. 

The above analysis and numerical results always used the 
peak amplitude of the ghost pulses and the signal pulses. 
The pulse width of the ghost pulses is ignored for simplicity. 
However, the IFWM ghost pulse broadens to \/3 times the 
signal pulse width [8], [9]. As the power is proportional to 
the pulse width, the ghosts pulse has VS times larger energy 
than the signal for the same peak amplitude. In the worst case, 
the X-axis of Fig. |3lmust scale by a factor of 3^/* = 1.32. 
However, the scale factor depends on the bandwidth of the 
optical and electrical filters in the receiver 

If optical match filter precedes the direct-detection DPSK 
receiver and the electric filter at the receiver has a wide band- 
width that does not distort the signal, the IFWM amplitude is 
increased by a factor of ■\/3/2 = 1.22 by the receiver, slightly 
less than the ratio of 1.32. If both the optical and electrical 
filters have a very wide bandwidth, allowing too much noise 
to the receiver, the peak amplitude directly transfers to the 
receiver. If the optical filter has a wide bandwidth but the 
electrical filter is a 0.75/T bandwidth Bessel filter, IFWM 
increases by a factor of 1.19 and 1.33 for To = 7.53 and 5 
ps, respectively. In practical system design. Figure |3] must be 



modified to take into account the design of both receiver and 
transmitter Note that the mean nonlinear phase shift of Fig. 
|3]is a simple system parameter to evaluate. 

Figure |3] shows that DPSK signal with IFWM can tolerate 
a far larger mean nonlinear phase shift of ($nl) than DPSK 
signal with nonlinear phase noise of ('I'nl) < 0.6 rad for 1- 
dB penalty [15]. However, Ho [15] is for return-to-zero (NRZ) 
signal without pulse distortion and deduces that RZ signal has 
lower tolerance to nonlinear phase noise. We are currently 
developing model for RZ signal with nonlinear phase noise in 
highly dispersive systems for a fair comparison. 

IV. Conclusion 

When the peak amplitude of IFWM induced ghost pulses is 
evaluated numerically, the error probability of DPSK signals 
can be found semi-analytically. For a SNR penalty less than 
1 dB, the mean nonlinear phase shift of the system must be 
less than 1 to 2 rad depending on the initial pulse width. An 
empirical approximation is also used to find the SNR penalty 
up to 2 dB. 
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